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ABSTRACT 

We investigate the possibility of correcting for the magnification due to grav- 
itational lensing of standard candle sources, such as Type la supernovae. Our 
method uses the observed properties of the foreground galaxies along the lines- 
of-sight to each source and the accuracy of the lensing correction depends on the 
quality and depth of these observations as well as the uncertainties in translating 
the observed luminosities to the matter distribution in the lensing galaxies. The 
current work is limited to cases where the matter density is dominated by the 
individual galaxy halos. However, it is straightforward to generalize the method 
to include also gravitational lensing from cluster scale halos. We show that the 
dispersion due to lensing for a standard candle source at z — 1.5 can be reduced 
from about 7 % to < 3 %, i.e. the magnification correction is useful in reducing 
the scatter in the Type la Hubble diagram, especially at high redshifts where 
the required long exposure times makes it hard to reach large statistics and the 
dispersion due to lensing becomes comparable to the intrinsic Type la scatter. 

Subject headings: gravitational lensing — supernovae: general 



1. INTRODUCTION 



Observations of high-redshift Type la supernovae (SNIa) have in the last decade or so 
lead to a dramatic paradigm shift in cosmology (Perlmutter et al. 1998; Riess et al. 1998; 
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Schmidt et al. 1998; Perlmutter et al. 1999; Knop et al. 2003; Tonry et al. 2003; Riess 
et al. 2004). Measurements of the luminosity distance to supernovae over a wide range of 
redshifts were used to break the degeneracy between cosmic fluids, as suggested by Goobar & 
Perlmutter (1995). The data clearly favors a universe dominated by repulsive dark energy, 
and which is presently undergoing accelerated expansion. The next step in observational 
cosmology is to test the nature of this dark energy, whether constant, i.e. compatible with 
Einstein's cosmo logical constant, or due to completely new physics. Observations of SNIa 
are among the leading astrophysical tools to explore this question further, as they probe 
the expansion history of the Universe directly. Large dedicated surveys are in progress (e.g. 
CFHTLS, ESSENCE, SDSSII) and even more ambitious space based projects are being 
planned for the future, e.g. the JDEM proposals, DESTINY, JEDI and SNAP. 

One thing in common for all these projects is the very large projected number of SNIa 
that eventually will populate the Hubble diagram used to derive cosmological parameters. 
Clearly, systematic uncertainties will (soon) become the limiting factor. While some of these 
uncertainties are due to our lack of knowledge of the SNIa physics and intrinsic properties, 
others stem from possible interactions of the supernova light (rest-frame UV and optical) 
near the source or along the line-of-sight (1-o-s), e.g. extinction by dust in the host galaxy or 
intergalactic medium. In this work, we focus on the gravitational interaction of photons along 
the 1-o-s, i.e. gravitational lensing. As supernova surveys become deeper, the measured source 
fluxes become increasingly more sensitive to the inhomogeneities in the matter distribution 
of the Universe. In Amanullah, Mortsell & Goobar (2003), the JDEM/SNAP mission was 
simulated using the SNOC Monte-Carlo package (Goobar et al. 2002) and it was found that a 
careful statistical treatment could be used to optimize the fitting of cosmological parameters 
from the Hubble diagram of SNIa taking into account the (asymmetric) redshift dependent 
lensing magnification distribution. Lensing on individual SNe have also been studied, in 
e.g Lewis & Ibata (2001); Mortsell, Gunnarsson k Goobar (2001); Bemtez et al. (2002); 
Gunnarsson (2004) by modeling the effect from the galaxies close to the 1-o-s to the SN. 

In this work, we investigate the accuracy to which lensing (de) magnification can be 
estimated on individual supernovae. For that purpose, we create mock galaxy catalogs 
with properties (e.g., galaxy magnitudes, redshifts and spectral types) based on luminosity 
functions derived from observations by Dahlen et al. (2005). Using the brightness of galaxies 
as a tracer of the gravitational fields along the 1-o-s, we use the multiple lens-plane package 
Q-LET (Gunnarsson 2004) to investigate the accuracy to which the magnification can be 
estimated as a function of the survey parameters, assumptions on M/L-ratios and halo 
shapes. In an accompanying paper (Jonsson et al. 2005), we apply the technique described 
here to investigate the lensing magnification probability distribution for 33 supernovae in 
the GOODS survey (Riess et al. 2004; Strolger et al. 2004). 
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In this paper, we will assume that dark matter halos in individual galaxy halos and 
small groups are most important for the lensing magnification of supernovae. For cluster 
size lenses, additional information is needed to model the gravitational potential, e.g., lensing 
of background galaxies. Though such a generalization of the method is straightforward, in 
the following we have assumed the uncertainty in the lensing magnification factor for the 
small fraction of SNe with foreground clusters to be of the same order as SNe with massive 
foreground galaxies. Also, we assume that the large scale dark matter structures in the 
Universe is traced by the luminous matter, i.e. that filaments and walls are populated by 
galaxies. This approach is reasonable as long as the luminous and dark matter are not 
anta-correlated, see §5. 

In §2, we discuss whether one needs to correct for lensing at all. Section §3 describes 
the underlying theory and §4 treats our method for estimating the accuracy of the lensing 
corrections. We summarize and discuss our results in §5. Throughout the paper we use 
natural units, where c = G = 1. We assume that the underlying cosmological parameters 
are Hq = 70 km s~ 1 Mpc _1 , the matter density £7m = 0.3 and the dark energy density 
f^A = 0.7. When no explicit redshift dependence is shown, quantities refer to present values 
(z — 0). Quoted magnitudes are in the Vega system. 



2. TO CORRECT OR NOT TO CORRECT 

Since the mean magnification due to gravitational lensing of a large number of sources 
is expected to be unity relative to a homogeneous universe, the question arises whether one 
should correct for the effect of gravitational lensing at all. 

Because flux, /, is conserved, it is the mean of the magnification factor, /z, that is 
equal to one, i.e. p, — 1, or defining fi = 1 + 5 where 5 is the fractional difference in 
luminosity from the unlensed (homogeneous universe) case, 5 = 0. The magnitude is given 
by m = —2.5 log / + const, and we can write 

m = m ° ~ hill) ln/i ' ^ 
where mo is the unlensed magnitude. Taylor expanding ln/i — ln(l + 8), we get 



2.5 

m = m °-hTTo 
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8 - 5 - + 0(8*) 



(2) 



with mean value m = m + 0.545 2 + 0(5 3 ). From this it is clear that the average lensed 
magnitude need not be equal to the unlensed magnitude. Note that in current surveys, this 
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effect is very small compared to, e.g., the intrinsic scatter of SN luminosities which is why 
the distinction is still unimportant. Note also that the mean magnification factor is unity 
only for random source positions. For an actual sample of observed SNIa, magnification bias 
can push the mean magnification to higher values. 

However, given that we have a sample of random source position SNIa and neglect the 
small corrections to m (or perform our cosmology fit using flux units), then m is an unbiased 
estimator for the population mean of the observed magnitudes 1 . Under these circumstances, 
neglecting the scatter due to lensing does not cause any bias in the fitted cosmological 
parameters and good statistics will help in beating down the error (e.g., Holz & Linder 
2004). There could still be good reasons to consider correcting for lensing effects. If we are 
able to reduce the scatter in the observed magnitudes and keep m as an unbiased estimator, 
then we are able to make more accurate cosmology fits. There are also cases where it is 
non-trivial to quantify the importance of the magnification bias, e.g., the case of SN 1997ff. 
In a similar context, the ability to correct individual lines-of-sight for gravitational lensing 
magnification would have a profound impact in our ability use gravitational wave "sirens" 
for measuring cosmological parameters, as their use as standard candles is ultimately limited 
by the lensing uncertainty (Holz & Hughes 2005). 

3. MODELING AN INHOMOGENEOUS UNIVERSE 

In this section we present a method to investigate the effects of gravitational lensing in 
an inhomogeneous universe. 

3.1. Halo Profiles 

Neglecting gravitational lensing is equivalent to assuming that matter is homogeneously 
distributed in the Universe. However, on small scales, the Universe is certainly inhomoge- 
neous. To investigate the effects of gravitational lensing of distant sources, a realistic model 
of the matter distribution in the Universe is needed. In the following, we describe how we 
(re) distribute the matter in our model universe using observations of the luminous matter. 

We assume that each galaxy is surrounded by a dark matter halo and that the mass of 
this halo can be estimated from the galaxy luminosity. However, inferring masses of dark 
matter halos from luminosities of galaxies is non-trivial. The effects of lensing by a halo 



As long as the variance of m is finite. 
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depends not only on its mass, but also on its density profile. Both the density profile and 
mass of dark matter halos are issues under debate. We have chosen to work mainly with two 
different halo models, Singular Isothermal Spheres (SIS) and the model of Navarro, Frenk 
and White (NFW; Navarro, Frenk & White 1997). 

The density profile of a SIS, psis(^) = o- 2 /(2irr 2 ), is characterized by its 1-o-s velocity 
dispersion a, which can be estimated from the galaxy luminosity via the Faber-Jackson 
(F-J) or Tully-Fisher (T-F) relations, approximately valid for elliptical and spiral galaxies 
respectively. Since the mass of a SIS halo diverges, msis(r) = 2a 2 r, we use a truncation 
radius r t . A commonly used scale for halo profiles in general is r 2 oo, defined as the radius 
inside which the mean mass density is 200 times the present critical density. For a SIS halo, 
r 2 oo and the corresponding mass within this radius, m 20 o, are given by 



\p2cj oto v^2i 



a 



3 
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The density profile of a NFW halo is 

r l = 

(r/r B )(l + r/r B ) 



PnfwM = 7r7: s - — A2 , (4) 



where r s is the scale radius for which approximately Pnfw oc r 2 and p s is the density at 
r ~ 0.5r s . 

The NFW halo is fully determined by m 20 o since r 20 o = [^2oo/(100i/o)] 1 ^ 3 an( i the scale 
radius r s and p s can be found numerically from m 2 oo (Navarro, Frenk & White 1997). 

In the following, we assume that the mass within r 200 is roughly the same for the SIS 
and NFW halo profiles, i.e. m 2 Jg = m^ Q W . We also set the truncation radius r t = r 20 o for 
both SIS and NFW halos. Varying r t does not alter gravitational lensing effects significantly, 
see §4.5. 



3.2. The Smoothness Parameter 

Very faint and/or small scale structures cannot all be seen in a magnitude limited survey. 
Also, for the method used in this paper, any matter not directly associated with individual 
galaxies such as completely dark halos and cluster halos need to be accounted for. In order 
to assure that the mean mass density in our model universe is kept constant, we keep the 
"remaining" mass, not accounted for when relating the dark matter to the luminous matter, 
as a homogeneous component. 

The homogeneous part can be characterized by the smoothness parameter r)(z), quan- 
tifying the fraction of smoothly distributed matter in our model universe (or our lack of 
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knowledge on the dark matter distribution in the real Universe). Since the fraction of galax- 
ies observed at a given magnitude limit is a function of redshift, and also since the Universe 
evolves, the smoothness parameter is expected to vary with redshift. 

The smoothness parameter in a given survey can be computed from the observed density 
of matter in clumps, i.e. in our case galaxies surrounded by dark matter halos, p s (z). If the 
redshift dependence of p g (z) can be factorized into a term (1 + z) 3 , scaling like the matter 
density, and an unknown factor f(z) originating from the magnitude limit of the survey and 
evolution, we can write 

p g (z)=p g (0)(l + z) 3 f(z). (5) 
Then the smoothness parameter is simply given by 

V(z) = 1 - ^f(z), (6) 

I'M 

where the density in galaxies at z = has been scaled with the present critical density to Qq. 
Once the galaxies have been associated with halos of definite masses, the comoving density 
of clumps as a function of redshift Qcf(z) can be estimated. We divide the distribution of 
galaxies into redshift bins and estimate QqJ(z) in each bin. The density of clumps in the 
i:th bin, centered on redshift Zi, is obtained through 

n G /(*) = (7) 

Pc * i 



V l = / A{ } , /9 AQdz, (8) 

Jz,-A Z /2 \n M a + z) 3 + n A } 1/2 



where rrij is the mass of a clump and p c is the critical density. The comoving volume of the 
i:th bin is given by 

r*+A*/2 D\{l + z) 2 

-az/2 [n M (i + z) 3 + n A ] 1/2 ' 

where Az is the width of the bin, AQ is the solid angle under study and Da is the angular 
diameter distance. Distances have been calculated using the angsiz routines described in 
Kayser et al. (1997), in which a smoothness parameter varying with redshift can be included. 
Note that the angular diameter distance Da used to determine the volume element in Eq. (8) 
above is calculated using the filled-beam approximation (77 = 1), since the volume is governed 
by the global expansion rate, which in turn is governed by the properties on very large scales 
where the Universe is homogeneous. 



3.3. Deriving Velocity Dispersions from Observed Luminosities 



Galaxy halo masses can be estimated from the velocity dispersion of galaxies. We 
calculate the velocity dispersion of each galaxy using absolute magnitudes (M B ) combined 
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with empirical F-J and T-F relations. For ellipticals, we use the F-J relation 

log 10 a = log 10 <t„ -—(M B - M* ) (9) 

7 

where M B is the characteristic magnitude and a* is the normalization in velocity dispersion. 
We use 7 = 4.4, as derived by Mitchell et al. (2005) using Sloan Digital Sky Survey (SDSS) 
data. To derive a*, we use equation (33) in Mitchell et al. (2005), where we use the relation 
M r = Mb — 1.32 to convert SDSS r-band magnitudes in AB system to standard S-band Vega 
normalized magnitudes. We have here assumed a typical color M B — M r = 1.20 for ellipticals 
in the AB-system, and an AB to Vega relation BAB=Bv e ga — 0.12. The normalization in 
velocity dispersion is given by 

l ogl0 a* = 2.2 - 0.091(M^ + 19.47 + 0.85z), (10) 

where we use M B = —21.04 derived for the early-type population by Dahlen et al. (2005). 
Equation (10) yields a* = 220 kms -1 at z — 0. Combining equation (9) and (10) gives an 
expression for the velocity dispersion 

log 10 a = -0.091(M B - 4.74 + 0.85z), (11) 

where we use z' — z for redshifts z < 1 and z' — 1 for z > 1. The redshift dependence of 
the relation accounts for the brightening of the stellar population with redshift. Since this 
evolution is poorly known at z > 1, we assume a flat evolution at these redshifts. As a 
measurement of the error in the derived relation, we use the observed scatter in the SDSS 
measurements by Sheth et al. (2003) 

rms(log 10 a) = 0.079[1 + 0.17(M B + 19.705 + 0.85z')], (12) 

where we again have transformed SDSS r-band to standard 5-band magnitudes. 

For the spiral and later type population, we use the T-F relation derived by Pierce & 
Tully (1992), with correction for redshift calculated by Bohm et al. (2004) 

log 10 V max = -0.134(M fl - AM B + 3.52), (13) 

where V max is the maximum rotation velocity for the galaxy. The correction due to redshift 
dependence is 

AM B = -1.22/ - 0.09. (14) 

The observed scatter in the relation derived by Pierce & Tully (1992) is rms(M B ) = 0.41, 
corresponding to 

rms(log 10 Knax) = 0.06. (15) 

At Mg, this is similar to the errors in the F-J relation above. Finally, the velocity dispersion 
in spiral galaxies is related to the circular velocity via o = V max / a/2. 
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3.4. Gravitational Lensing with Multiple Lenses 



A typical source 1-o-s within some angular radius 8 S will contain more than one lens. This 
requires the multiple lens-plane algorithm (see Schneider, Ehlers & Falco 1992; Gunnarsson 
2004, for further details), which takes into account each lens along the 1-o-s by projecting 
the lens' mass distribution onto a plane and then traces the light-ray from the image plane 
(first lens-plane) back through all lens-planes up to the source plane where the magnification 
and intrinsic position can be found. 



In the following we denote angular diameter distances between redshifts z { and Zj by 
Dij. We use o for observer, s for source and d for lens (deflector). When z% = 0, that index 
is omitted. 



Each halo is truncated in 3D at r = r t , then, upon projection onto a plane, the corre- 
sponding surface mass density will smoothly go to zero at the projected truncation radius. 
The projection can be done analytically for our lens models. For simplicity, we start by 
considering a single lens-plane. The equations can be simplified if we let £ be the impact 
parameter on a halo and define x = £/£o and x t = r t /£o, where 



for the SIS and 



4na 2 D d D ds 



£o = r B 



for the NFW halo. Then, the projected density n(x) can be written as 
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x < 1 < x t 

1 < X < x t 

(20) 
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and 



. . 2 / x^ "I - 2 i 

«NFW(1) = TT^s ~ ■ (21) 



3 v 

Here, k s = p s r s /£ cr , where S cr = D s / AnD&D&s is a critical density related to strong lensing. 
Note that x t > 1 must be assumed for the NFW and that k = for x > x t for both halo 
types. 

The general expression for the deflection angle for circularly symmetric lenses is 

2 f x 

a(x) = sa(x) = s— / x ' K(x')dx ', (22) 
x Jo 

where s = £,oD s /D d D ds , resulting in 

„(*> = I (^EZ) + ii^rzzj m 

for the SIS model. For the NFW halo, numerical evaluation is needed. 

As the magnification factor, //, for all halo models is obtained with distances calculated 
with the ^-dependent rj function, this is the universe relative to which // is found (implying 
// > 1 for primary images). In the following, we will quote magnifications, fi, relative to 
a universe with homogeneously distributed matter, the filled-beam value (fb) where p, — 1. 
The magnifications are related by 



Ifb x 2 



4. SIMULATED SURVEYS 

In order to study gravitational lensing corrections, we perform Monte Carlo simulations 
where we calculate the magnification factor for random source positions in mock galaxy cat- 
alogs. By varying the assumptions of the galaxy mass distributions as well as the magnitude 
limit of the observations, we can estimate the accuracy to which it is possible to correct for 
the lensing magnification. 

For all lensing calculations we have used the publicly available f ortran 77 code Q-LET 2 
(Gunnarsson 2004), although substantially modified. The code fully utilizes the multiple 



2 Available at http://www.physto.se/~cg/qlet/qlet.htm 
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lens-plane algorithm and has been used previously by Gunnarsson (2004) and Riess et al. 
(2004) to study lensing effects on supernovae. 

As our simulation base, we create for each Monte Carlo realization a mock galaxy 
catalog designed to reflect the distribution of galaxies expected in a circular cone around a 
random 1-o-s. To characterize the galaxy population we use the 5-band rest-frame Schechter 
luminosity function (LF) derived by Dahlen et al. (2005) using GOODS CDF-S observations. 
The LF is used to generate the number of expected galaxies within the cone where we take 
into account Poissonian fluctuations but do not include effects of galaxy correlations or 
cosmic variance. The same LF is used to assign absolute magnitudes to each object within 
the range —23 < M B < —16. To account for evolutionary effects, we include a brightening 
of the 5-band characteristic magnitude by ~1 mag to redshift z = 1 as discussed in §3.3. 

A random spectral type is assigned according to the type-specific LF of early-types, 
late-types and starburst galaxies at z ~0.4 in Dahlen et al. (2005). We thereafter assign 
early-type galaxies an elliptical morphology and late-type galaxies a spiral morphology. We 
assume that the fraction of galaxies with elliptical morphology is constant over the redshift 
range investigated. 

The redshift of each object is assigned with a probability proportional to the volume 
element, dV(z)/dz, which is equivalent to assuming a constant comoving number density of 
galaxies with redshift, i.e. we do not include any evolution of the number densities due to 
e.g., mergers or large scale structures. The galaxy is finally given a random position within 
the 1-o-s cone. 

Besides redshift, absolute magnitude, spectral type and position, we also calculate the 
apparent magnitude in the observed /-band for each object. This allows us to draw subsam- 
ples from the catalog with imposed magnitude cutoffs as is the case for real observations. 
Furthermore, to resemble an observational situation where redshifts are determined photo- 
metrically, we also have the option to add a random error to the redshifts. These errors are 
calculated using simulations and depend on redshift, spectral type and detection S/N. For 
bright objects with high S/N, errors are typically A z = (|-2 phot — 2 true |/(l + ^spcc)) ~ 0.05, 
while errors for faint objects (mostly at high z) can be as large as Az ~ 0.3. The errors for 
early-type galaxies are typically a factor two smaller compared to late-type galaxies when 
comparing at the same apparent magnitudes. Besides this "Gaussian" contribution to the 
error distribution, a fraction of the galaxies may also be assigned "catastrophic redshifts" 
with large errors. We discuss this further in §4.4. 

Figure 1 shows the simulated accuracy of the photometric redshifts for a survey with 
limiting magnitude I < 27 (S/N=10). The bottom panel shows the generated (input) 



redshifts vs. the photometric redshifts, while the top panel shows the difference between 
generated and photometric redshifts as a function of galaxy magnitude. 

4.1. Understanding the Magnification Uncertainties 

We have identified and addressed the following uncertainties when estimating the lensing 
magnification of a specific source given a galaxy catalog: 

• Finite field size 

• The intrinsic scatter in, and accuracy of, the F-J and T-F relations 

• Redshift and position uncertainties 

• Choice of halo profile 

• The magnitude limit 

These sources of error are addressed individually in the following sections. Our reference 
model consists of NFW halos truncated at r t = r 2 oo, velocity dispersion/circular velocity 
normalizations of 220 kms -1 for ellipticals, 203 kms" 1 for spirals and source redshift z = 1.5. 
In Figure 2, the Probability Distribution Function (PDF) for lensing magnifications for the 
reference model is shown. The lensing dispersion at z — 1.5 is ~ 7%. We analyze the results 
by comparing the distribution of magnifications in the reference model with the distribution 
obtained after performing a correction with the above-mentioned uncertainties. We denote 
the uncorrected value /x re f and the corrected one g M , where 

/^ref 



where // est is the estimated magnification factor including one or more uncertainties. Cor- 
rections will reduce the uncertainties from lensing whenever the width of the distribution 
of is smaller than the corresponding width in the /i re f distribution. If no uncertainties 
were present /i est = // re f and — 1 implying a perfect correction. As a measure of the 
width we give the standard deviation of the distribution. Since many of the distributions 
are non-Gaussian, we also report the 68 % and 95 % confidence levels. 



(25) 



4.2. Finite Field Size 



The mean magnification relative to an homogeneous universe of a large number of 
sources lensed by randomly distributed matter is expected to be unity due to photon number 
conservation. When we model a lensing system, only galaxies within angular radius 6 S of the 
position of the source on the sky are taken into account and thus p, < 1. If 9 S is increased, 
more lenses are added and the mean magnification increases. This dependence is illustrated 
in Figure 3, showing the mean magnification of 5000 point sources, as a function of 9 S . 
The mean magnification increases rapidly for small 6 S , but only slowly for 6 S larger than an 
arc-minute, where the error is ~ 1 %. In our simulations, we use 60" as a cutoff to save 
computing time. In a real survey, a cutoff will have to be introduced for practical purposes 
since only a limited portion of the sky will be observed. In order to avoid a systematic 
bias due to the finite field size for a given survey, the computed magnifications should be 
corrected with a factor corresponding to the inverse of the mean magnification for the cutoff 
radius used. In the following, we have neglected this small (~ 1 %) correction. Furthermore, 
going to larger 6 S would not render fx being exactly unity since some flux is lost whenever 
multiple imaging occurs. Q-LET gives the magnification and intrinsic source position of a 
given image, not observed position and magnification of a given source. Therefore, in the 
rare cases of multiple imaging, only one of the images will be taken into account resulting in 
some flux loss. Note also that since random 1-o-s and random source positions are different 
(e.g., Schneider, Ehlers & Falco 1992), we have to use a magnification dependent weighting 
procedure to see whether each simulated event should be kept or discarded in order to get a 
sample of random source positions (see e.g., Goobar et al. 2002). 



4.3. The F-J and T-F Relations 

The F-J and T-F relations give the velocity dispersions of the luminous matter and 
we make the assumption that the dark matter that constitutes the halo follows scales in the 
same way. Both the Faber- Jackson and the Tully-Fisher relations have an intrinsic scatter 
with an rms estimated in §3.3. To study the effect of this scatter, we add random offsets to 
the F-J and T-F relations when calculating the halo mass. Since we do not want to bias 
the total mass in our simulations, we distribute the offsets using a Gaussian distribution in 
a 3 (since mass oc a 3 ). We derive the width of the Gaussian (one sigma value in a 3 ) from the 
rms in logi (<r) and logi (V r max ) using Eqs. (11)-(12) (F-J) and Eqs. (13)-(15) (T-F). 

In panel a) in Figure 4, we compare the distribution of the corrected value due to the 
scatter in the F-J and T-F relations with the distribution of magnifications in the reference 
model (dashed line). Note that if we knew all velocity dispersions and circular velocities 
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exactly, would be represented by a 5-function at — 1. We see that the intrinsic scatter 
in the velocity dispersion and circular velocity causes a dispersion of g M of approximately 
3%, a factor of ~ 2.6 less than the original dispersion. Panel a) in Figure 5 shows /i re f — 1 
vs g M — 1 for each individual source. For 82 % of the sources, the corrected luminosity will 
be better than the uncorrected one. 

Besides the intrinsic scatter in the F-J and T-F relations, there is also a possible 
systematic uncertainty in a*. To estimate this, we use SDSS data in Bernardi et al. (2003) 
where photometric and spectroscopic parameters for a sample of ~ 9000 early-type galaxies 
in the redshift range 0.01 < z < 0.3 are given, including K-corrections and accurately 
measured velocity dispersions. We fit a straight line Mb = b log a + a and estimate the error 
in a* for a given M B by propagating the errors in the parameters a and b. For M* B = —21.04, 
we obtain a* = 218 ± 7 kms -1 . In panel b) and c) in Figure 4, we investigate the effect 
of a systematic shift of ±10 kms -1 in er*. The dispersion in due to such a shift is quite 
small or at the order of 1 %. Panels b) and c) in Figure 5 shows /z rc f — 1 vs — 1 for each 
individual source. The corrected value will be better than the uncorrected for > 95 % of the 
sources. 

A further source that may increase the scatter in the magnification is the possible 
misclassification of galaxy morphology. An elliptical galaxy wrongly classified as a spiral, 
leads to an underestimation of the underlying mass, and vice versa if a spiral is misclassified 
as an elliptical. To investigate the possible effect of this, we first use a set of simulated 
galaxies with known morphological types (i.e. an exponential radial profile for spirals and 
a de Vaucouleurs profile for ellipticals) and measure how many are correctly recovered in 
an observational setup resembling the GOODS. To classify galaxies, we use the GALFIT 
software (Peng et al. 2002), which measures the slope of the radial profile and therefore 
allows a discrimination between ellipticals and spirals. At a S/N=10 detection limit (m~25), 
we find that ~25 % of the ellipticals and ~8 % of the spirals are misclassified. The fraction 
of misclassified galaxies quickly drops to ~1 % at a magnitude ~2 mag brighter than the 
detection limit. We then use these results to estimate the effect of misclassification on the 
derived magnifications. We find that the increase in the scatter in the magnification is 
~0.5 % due to this effect. Therefore, misidentification of galaxy morphology should only 
affect the results marginally. However, for ground-based surveys with low resolution, the 
effect may be larger. 
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4.4. Redshift and Position Uncertainties 

Uncertainties in the redshifts of the lenses will alter the results both through the un- 
certain distances between different lens-planes and by introducing an uncertainty in their 
absolute magnitudes used in the F-J and T-F relations. 

In an ideal observational situation, all redshifts are determined spectroscopically. In 
many real situations, however, only photometric redshift are available due to, e.g., the faint- 
ness of the galaxies and the large number of sources. To investigate the effects of photometric 
redshift uncertainties, we add a random offset to the redshift of each object. The size of the 
offset depends on redshift, apparent magnitude and spectral type of the object and is drawn 
from the simulated error distribution discussed above. 

In panel d) in Figure 4, we show the distribution of the corrected value due to a 
random offset to the redshift of each lensing object. The induced error in the estimated 
magnification is less than 1 %. The corresponding panel in Figure 5 shows /i rc f — 1 vs g M — 1 
for each individual source. In this case, 96 % of the sources have corrected values which are 
better than when uncorrected. 

Besides the Gaussian-like distribution of the photometric redshift errors investigated 
above, there is also a possibility that a fraction of the objects get "catastrophic redshifts" 
with large errors, so called outliers. E.g., by comparing with ~1400 spectroscopic redshifts 
in the CDF-S and HDF-N, we find that the GOODS photometric redshifts have about 3% 
outliers with A z > 0.3. The redshift probability distribution, for a majority of these objects 
are characterized by a primary (~Gaussian) peak combined with a less pronounced secondary 
peak. Outliers are foremost objects assigned the redshift of the primary peak, but where 
the true redshift is that of the secondary peak. To estimate the effect of outliers, we use the 
galaxies in the GOODS and simulate the case where we distribute the photometric redshifts 
over the full probability distribution, including the secondary peak. This will allow ~3 % 
outliers. We compare this with the case where we only include redshifts in the primary peak 
(i.e. objects with A z < 0.3). We find that the increase in lensing dispersion due to the 
population of outliers is less than 1 %. One reason for the small effect is that outliers are 
mainly faint and very blue objects, which therefore should have relatively small masses. 

Any error in the exact positions of the lensing galaxies can also affect the resulting mag- 
nification. Apart from the observational error, such an effect can be due to a misalignment 
between the luminous and the dark matter in a given galaxy. We have investigated the effect 
of a Gaussian random shift with o"p OS = 0.5 arcseconds of all lensing galaxies along the 1-o-s. 
Even such a large shift of all galaxy positions results in a distribution of corrected values 
with a dispersion of less than 0.5 %. 
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4.5. Choice of Halo Profile 

The choice of halo model is only important in those lens-planes where the light-ray 
passes through a halo. If passing outside, the halo will act as a point mass and when 
m tot° = m loo = m 2oo W the two different halo models give exactly the same results. 

We have performed simulations where all halos were of SIS instead of NFW type. The 
effect of different halo profiles are also present in the realistic and pessimistic case simulations 
below. 

Panel e) in Figure 4 shows the PDF of g M when assuming SIS halos instead of NFW 
as in the reference model. The dispersion is less than 1.5%. In > 90 % of the cases, the 
corrected value is better than the uncorrected, see panel e) in Figure 5. 

We have also investigated how important the assumption on the truncation radius is for 
the resulting magnification distribution by running tests with 0.75 x r 20 o < r t < 1.25 x r 20 o 
for the SIS model. The uncertainty in the resulting magnification gives a distribution of 
corrected values with a dispersion of ~ 0.5 — 1 %. Since the NFW profile falls of as 
Pnfw oc r^ 3 at large radii as compared to psis oc r -2 for SIS halos, this should be considered 
a very conservative limit on the effect of changing the truncation radius. 



4.6. Magnitude Limits 

Our reference model uses a constant comoving mass density of galaxies, implying a con- 
stant smoothness parameter, 77. In a real scenario with an observational magnitude limit, an 
increasing fraction of galaxies drop out at higher redshift. The faint high redshift galaxies 
will not be seen and hence not included as lenses in the magnification calculation but instead 
added as homogeneously distributed matter. Therefore, when deriving the smoothness pa- 
rameter from observations, r)(z) increases with redshift even if the 'underlying' smoothness 
parameter is constant. For each simulation with a finite magnitude limit, a new 77-function 
is computed using the method described in §3.2. 

For I = 27 mag the distribution of is very narrow and for I = 29 mag, it is in 
principle a 5-function. Panel f) in Figure 4 shows the PDF of q^ for I = 23. The dispersion 
is ~ 2 %. In 86 % of the cases, the corrected value is better than the uncorrected, see panel 
f) in Figure 5. 

In Figure 6, we compare q^ — 1 as a function of source redshift for models with magnitude 
limits / = 23 and / = 25 with /j re f — 1. Even for source redshifts as high as z ~ 2, a magnitude 
limit of I = 23 does not significantly impair our results. 
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Photometric errors in the apparent magnitudes translates to an increased scatter in the 
absolute magnitudes and therefore also in the derived velocity dispersions and masses. At 
the faintest magnitude limits considered here, S/N=10, typical errors are ~ 0.1 mag. For 
ellipticals, this corresponds to an increased dispersion of AlogioO" ~ 0.01 (using Eq. 10). 
This is significantly less than the intrinsic scatter in the F-J relation which is rms(logio)o" ~ 
0.08 (at ~ M|j, Eq. 12). For the T-F relation the increased scatter due to photometric 
errors is AlogioVmax ~ 0.013 (using Eq. 13), again significantly less than the intrinsic scatter 
rmsQogioKnax) ~ 0.06. So even for the faintest galaxies considered, the errors in apparent 
magnitude should not affect results more than marginally. 



4.7. Realistic and Pessimistic Scenarios 

We have studied the uncertainty in the lensing correction in a realistic scenario where a 
reasonable error budget is assumed. In this case we have assumed 50% NFW, 50% SIS, no 
central value shift of the velocity dispersion normalization but a dispersion around this value 
and a magnitude limit of / = 25. Lens redshifts were assumed to be distributed around their 
reference values. As the correct model we have used the reference NFW model as above. 

A pessimistic scenario for space based surveys has also been studied where we have 
maximized the uncertainties in relating the luminous and the dark matter. However, one 
could easily imagine worse cases in a ground-based experiment. For our scenario the erro- 
neous assumptions were: SIS halos, central value of velocity dispersion normalization shifted 
+ 10 kms -1 and distributed around this value, a magnitude limit of / = 25, lens redshifts 
assumed to be distributed around their reference model values, an offset in lens positions 
as described in §4.4 and finally a truncation radius of 1.25 x r 2 oo- We consider this being a 
pessimistic but not completely unrealistic scenario. 

The left panels in Figure 7 show results for the realistic scenario, the right panels for 
the pessimistic scenario. For the realistic case, has a dispersion of ~ 3 % and (q^ — 1) < 
(/W — 1) f° r ~ 80 % of the sources. In the pessimistic case, the corresponding numbers are 
~ 3% and ~ 77%, i.e. our ability to correct for lensing is more or less unimpaired when 
going from a realistic to a pessimistic scenario. The bottom row shows q^ — 1 as a function 
of source redshift for the two scenarios. The confidence levels of the realistic case vs z can 
be well fitted with straight lines and these are found in Table 1 expressed in magnitudes. 
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5. SUMMARY AND DISCUSSION 

We have investigated the accuracy to which lensing magnification can be estimated on 
individual lines of sight using the observed properties of the foreground galaxies of each 
source. 

The result depends on the uncertainties in translating observed galaxy luminosities to 
the (invisible) matter distribution in the lensing galaxies. We have shown that none of the 
studied uncertainties neither individually nor combined will render the corrected distribution 
of magnifications wider than the dispersion from lensing. Even for a pessimistic scenario, 
the dispersion due to lensing for a standard candle source at z = 1.5 can be reduced by 
a factor > 2, comparable to the result for a realistic scenario. The reason our pessimistic 
case result is not significantly worse than for the realistic case, is that the uncertainties are 
dominated by the scatter in the F-J and T-F relations for both scenarios 3 . At lower redshifts 
(z < 0.5), the effects from lensing are small and correcting for lensing is not likely to improve 
the results. 

Even though the fraction of SNe lines-of-sight passing through galaxy cluster lenses is 
expected to be relatively low, these are potentially important due to magnification bias and in 
surveys specifically aimed at using cluster potentials as gravitational telescopes (Gunnarsson 
& Goobar 2003). In those cases, individual modeling of the cluster potentials using, e.g., weak 
and strong lensing of background galaxies is needed to correct the observed magnitudes for 
the lensing magnification. Alternatively, one can choose to discard SNe background to galaxy 
clusters with very uncertain matter distributions when estimating cosmological parameters. 

Our method also takes into account gravitational lensing from large scale dark matter 
structures such as filaments and walls as long as the matter density is dominated by the 
individual galaxy size halos. If we very conservatively assume that large scale structures 
are completely uncorrelated with the luminous matter, we expect the lensing magnification 
contribution to be less than 2 % on scales larger than 5 arcminutes (for a source redshift of 
unity) (Cooray et al. 2005). 

For a given galaxy catalog, the computed magnifications do not depend strongly on 
the cosmological parameters used. However, since the formation of matter structure is a 
function of cosmology, it should in principle be possible to determine the cosmology from the 
observed distribution of standard candle luminosities. In this case, the use of magnifications 



3 We were able to increase the width in the distribution for the pessimistic case scenario with ~ 50 % 
by modelling 20 % of the galaxies as point masses. However, we consider such a compact mass distribution 
at galaxy scales too contrived to be included in the simulations. 
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of e.g. SNIa would be similar to using shear of background galaxies as in weak lensing studies. 
Such a study would require large statistics of very well observed SNIa and will probably have 
to await future dedicated missions such as the proposed SNAP satellite 4 . 

For current high-z SNIa observations the concern is to correct for the magnification 
and investigate the possibility for magnification bias. Such a study for SNe in the GOODS 
fields is described in a accompanying paper (Jonsson et al. 2005) where magnification bias is 
shown to be negligible but lensing for individual SNe can be estimated quite robustly from 
foreground galaxy observations and thus be corrected for. In fact, as long as the luminous 
and dark matter are not anti-correlated, we would expect to be able to reduce the scatter 
in the Hubble diagram by assuming that dark matter follows light. Thus, we find that even 
though the exact relation between luminous and dark matter is uncertain, correcting for 
gravitational lensing using observed galaxy properties should be harmless at the worst and 
very useful at the best. 
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Equation: y = kz + m 



68 % magn. 
68% demagn. 
95 % magn. 
95% demagn. 



Uncorrected 
k m 
-0.038 0.005 

0.060 -0.009 
-0.140 -0.003 

0.086 -0.017 



Corrected 
k m 
-0.020 0.008 
0.016 0.006 
-0.055 0.008 
0.035 0.012 



Table 1: Fits to the confidence levels of Fig. 7 expressed in magnitudes. Note that these are 
only valid for z > 0.25. However, below this value the corrections are indeed very small. 
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Generated z 



Fig. 1. — Bottom panel shows generated (input) redshifts vs. derived photometric redshifts 
for a simulated survey with limiting magnitude I < 27 (S/N=10). Top panel shows the 
difference between generated and photometric redshifts as a function of observed galaxy 
magnitude. 
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Fig. 2. — PDF of magnifications for the reference model with NFW halos truncated at 
f\ = ^200) no magnitude limit, velocity dispersion/circular velocity normalizations 220 kms -1 
for ellipticals, 203 km s" 1 for spirals, no uncertainty in halo redshifts and source redshift 
z = 1.5 




Fig. 3. — Mean magnification /x of 5000 lines-of-sight as a function of 9 S . The sources all 
have redshifts z = 1.5. 
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Fig. 4. — PDFs of where the following uncertainties have been studied: a) The scatter in 
the velocity dispersion of the F-J and T-F relations, b) Offset in central value of <r* of +10 
kms -1 , c) Offset in central value of a* of —10 kms -1 , d) Lens redshifts distributed around 
their original values, e) Assuming SIS halos instead of NFW and f) Assuming a magnitude 
limit of / = 23 mag. The PDF before correction (dashed line) is also plotted for reference. 
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Fig. 5. — Spread in /x re f — 1 vs spread in — 1 where the following uncertainties have been 
studied: a) The scatter in the velocity dispersion of the F-J and T-F relations, b) Offset in 
central value of cr* of +10 kins" 1 , c) Offset in central value of a* of —10 kins" 1 , d) Lens 
redshifts distributed around their original values, e) Assuming SIS halos instead of NFW 
and f) Assuming a magnitude limit of / = 23 mag. The gray shaded areas show where the 
value after correction is closer to reality than before correction. The figures in the upper 
left corners give the percentage of events residing in the gray areas, where the correction 
improves the resulting accuracy of the derived distance. 
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Fig. 6. — — 1 as a function of source redshift for a models with magnitude limits / = 23 
and I = 25, using 1700 lines-of-sight for each redshift. The gray shaded areas show the 
68% and 95% confidence levels of // re f — 1, the distance from the correct value 1 before 
correction. The dashed and dotted lines show the same confidence levels after correction 
clearly indicating the improvement. 
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Fig. 7. — Left column shows a realistic scenario with 50 % NFW, 50 % SIS, values of velocity 
dispersion distributed around their central value, lens redshift values distributed around their 
reference model values and finally using a magnitude limit of / = 25 mag. The right column 
shows a pessimistic scenario with with 100 % SIS, values of velocity dispersion normalization 
distributed around a +10 kms -1 shifted central value, a small random offset (<7 pos = 0.5") 
in lens positions, truncation radius of 1.25 x r 2 oo and the rest as in the realistic case. The 
upper row gives the PDFs and confidence levels for sources at redshift z = 1.5, middle row 
shows the spread in /x ref — 1 vs the spread in — 1. The bottom row shows confidence levels 
of the two scenarios vs source redshift. 



